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Abstract. We compute the order corrections to the one particle inclusive electroproduction cross 
section of hadrons with non vanishing transverse momentum. We compare our results with HI data 
on forward production of 7t°, and conclude that the data is well described by the DGLAP approach, 
within the theoretical uncertainties. 



INTRODUCTION 

The precise measurement of final state hadrons in lepton nucleon deep inelastic scat- 
tering constitutes an excellent benchmark for different features of perturbative quantum 
chromodynamics. Among them, the calculation of higher order corrections, which have 
been explored and validated for most processes up to next to leading order (NLO) accu- 
racy. For the one particle inclusive processes only very recently there has been progress 
beyond the leading order (LO) JE & S 0> and up to now there were no analytic 
computation of the &(ocf) corrections for the electroproduction of hadrons with non 
vanishing transverse momentum. The analytic computation of the 0(oc}) corrections al- 
lows us to check factorization in a direct way, which means that collinear singularities 
showing up in the partonic cross section factorize into parton densities (PDFs) and frag- 
mentation functions (FFs). As a consequence of this explicit cancellation, the resulting 
cross section is finite and can be straightforwardly convoluted with PDFs and FFs in a 
fast and stable numerical codes. The analytical result is still sufficiently exclusive and 
keeps the dependence on the rapidity and the transverse momentum of the produced 
hadron, allowing a detailed comparison with the experimental data. In the following we 
summarize the results obtained in Ref. [6] 



^(af ) QCD CORRECTIONS 

We consider the process 

l(l)+P{P)^l'(l')+h{P h )+X, (1) 

where a lepton of momentum I scatters off a nucleon of momentum P with a lepton 
of momentum /' and a hadron h of momentum Pj r tagged in the final state. Omitting 



target fragmentation at zero transverse momentum, which has been discussed at length 
in QlJ, |2|] , the cross section for this process can be written as 
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where <xX is the partonic level cross section corresponding to the process and is cal- 
culated order by order in perturbation theory through the related parton-photon squared 
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matrix elements (i, j) for the i + y — > j + X processes 
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in terms of the standard kinematical variables @]. 

At order-a A 2 , the partonic cross sections receive contributions from the following 
reactions: 

Y+q{q) 
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Real contributions 



Virtual contributions 
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where any of the outgoing partons can fragment into the final state hadron h. 

At variance with the \pr \ = case, where the integration over final states leads to 
overlapping singularities along various curves in the residual phase space, here the only 
remaining singularities are found at z = and thus they can be dealt with the standard 
method. After combining real and virtual contributions to a given partonic process, the 
cross section can be written as 
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+ Clf(p,y,z) + 0(e) 



(5) 



(21 

where the coefficient of the single poles, ^j ; j (p,j,z), as well as the finite contributions 

(2) 

C\j (p ,y,z), include 'delta' and 'plus' distributions in z. The IR double poles present 
in the individual real and virtual contributions cancel out in the sum, providing the 
first straightforward check on the angular integration of real amplitudes and the loop 
integrals in the virtual case. In the real terms, the above mentioned double poles come 
from the product of a pole arising in the integration over the spectators phase space and 
a single pole coming from the expansion of z~ 1+e factors. Double poles in the virtual 
contributions always arise from loop integrals. 

The remaining singularities, contributing to the single pole, are of UV and collinear 
origin. The former are removed by means of coupling constant renormalization, whereas 
the latter have to be factorized in the redefinition of parton densities and fragmentation 
functions. 



PHENOMENOLOGY 



In Figure [T] we show the LO and NLO predictions for the electroproduction of neutral 
pions as a function of xb and pj, respectively, in the kinematical range of the HI 
experiment [J, together with the most recent data for the range pj > 3.5 GeV. The 
cross sections are computed as described in the previous sections, applying HI cuts 
and using MRST02 parton densities [8]. Similar results are found using other sets of 
modern PDFs. For the input fragmentation functions, we use two different sets, the ones 
from reference ^ denoted as KKP and those from Jioll referenced as K. We set the 
renormalization and factorization scales as the average between Q 2 and p\, and we 
compute a s at NLO(LO) fixing Aqcd as in the MRST analysis. The plots clearly show 
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FIGURE 1. LO and NLO cross sections, including experimental cuts as explained in the text, as a 
function of xg. HI data |7] for the range pj > 3.5 GeV are also shown. 

that the NLO cross sections are much larger than the LO ones, even by the required 
order of magnitude in certain kinematical regions. Another interesting feature is that the 
uncertainty due to the choice of a fragmentation functions set is also quite noticeable, 
this fact driven by the different gluon content of the two sets considered here. Low Q 2 
bins seem to prefer KKP set, which have a larger gluon-fragmentation content, whereas 
for larger Q 2 both sets agree with the data within errors. LO estimates show a much 
smaller sensitivity on the choice of fragmentation functions, since gluon fragmentation 
does not contribute significantly to the cross section at this order. 

The rather large size of the K-factor can, then, be understood as a consequence of the 
opening of a new dominant ('leading-order') channel, and not to the 'genuine' increase 
in the partonic cross section that might otherwise threaten perturbative stability. The 
dominance of the new channel is due to the size of the gluon distribution at small xb and 
to the fact that the HI selection cuts highlight the kinematical region dominated by the 
Y+g — ► g + q + q partonic process. 

In Figure El we show the different contributions to the cross section discriminated by 
the underlying partonic process. Notice that at very small xb the gg term can be by itself 
several times larger than the LO contribution, remaining larger or comparable even for 





higher xb values. The forward selection is also responsible of the scale sensitivity of the 
cross section, as it supresses large components with small scale dependence whereas it 
stresses components as gg whose scale dependence would be partly canceled only at 
NNLO 
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